Gating of a molecular transistor: Electrostatic and Conformational 
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We derive a general result that can be used to evaluate and compare the transconductance of dif- 
ferent field-effect mechanisms in molecular transistors, both electrostatic and conformational. The 
electrostatic component leads to the well-known thermal limit in the absence of tunneling. We show 
that in a standard three-terminal geometry and in the absence of strong electron-phonon coupling, 
the conformational component can lead to significant advantages only if the molecular dipole mo- 
ment ^ is comparable to etox, tox being the thickness of the oxide. Surprisingly this conclusion is 
independent of the "softness" of the conformational modes involved, or other geometrical factors. 
Detailed numerical results for specific examples are presented in support of the analytical results. 

PACS numbers: PACS numbers: 85.65.+h, 73.23.-b,31.15.Ar 
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Device miniaturization is progressively heading to- 
wards solid state electronic components that are molecu- 
lar in nature Molecular transistors are of great cur- 
rent interest from both basic and applied points of view. 
Despite theoretical proposals H, HI and experimental 
reports 0,0,13 of three-terminal molecular devices, their 
general physical principles are not yet well understood. 
In a standard silicon MOSFET the gate modulates the 
current by controlling the channel charge through its elec- 
trostatic potential. Good transistor action requires that 
this channel potential respond much more to the gate 
than to the drain, implying an oxide thickness iox that 
is much smaller than the channel length L (Fig. ^). A 
molecular transistor operating on electrostatic principles 
ought to have the same design limitation, so that even 
nominal gate control in a 10 A molecule such as phenyl 
dithiol (PDT) 0| demands an oxide that is prohibitively 
thin Q. It is therefore natural to investigate alternate 
principles of molecular transistor action, such as by utiliz- 
ing conformational degrees of freedom 4]. For instance, 
a gate field coupled with a molecular dipole moment fl 
could cut off the current by tilting it away from a contact 
(Fig. nj)), or by twisting one part relative to the other, 
breaking its conjugation [iflllllj (Fig.^^)- Such a mech- 
anism could respond more strongly to the gate than the 
drain if fl is engineered to lie along the source-drain field. 

In this paper, we derive a general expression 

for the transconductance per unit current Qm/I = 

{l/I){dI/dVg) {I: current, Vg-. gate voltage) for dif- 
ferent transistor mechanisms, both electrostatic and con- 
formational ^3 . For a large separation of time scales be- 
tween electronic and vibronic modes, as is typical in room 
temperature molecular conductors [igl, the electrostatic 
and conformational transconductances are additive, with 
a well-known electrostatic maximum: 



e/3 = {25mVy 



(1) 



{j3: inverse thermal energy), achieved only if tox ^ L and 
if there is negligible tailing in the density of states. In ad- 
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FIG. 1: In a standard transistor (a) the gate voltage Vg con- 
trols the current / by controlling the induced charge. But 
in conformational transistors the gate controls the current 
through the configurational parameter 9. This could be ei- 
ther (b) through a variation in coupling with the contact, or 
(c) a variation in transmission through the molecule itself. 



dition, we find a fundamental conformational maximum, 
independent of modal stiffness or molecular geometry: 



[9T'/l]^,^^ePit,/et,^). 
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These simple expressions for gate control are later 
corroborated with detailed numerical results that self- 
consistently combine a semi-empirical description of the 
device Hamiltonian with a nonequilibrium Green's func- 
tion (NEGF) description of transport. Based on these 
results, it seems that good molecular transistor action in 
a standard MOS geometry with a modest gate insulator 
thickness would require degenerately doped semiconduct- 
ing contacts, along with the incorporation of a molecular 
dipole dipole aligned along the source-drain direction and 
large enough to overcome any room-temperature thermal 
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averaging over conformations, ie, larger than etox- 

We assume above that conducting and non-conducting 
molecular states lie in the same valley in the conforma- 
tional potential landscape. The gate voltage shifts the 
minimum within this valley but thermal excitations lead 
to a statistical average over the entire valley, making it 
difficult to change the current any faster than the fun- 
damental limit implied by Eq. |21 One could get around 
this limit if the conducting and non-conducting molecular 
states belong to distinct metastable valleys |j| , such that 
thermal averaging takes place only over one valley or the 
other. The role of the gate then is to shift the molecule 
from one valley to the other through a large impulsive 
force. Such a mechanism, however, is very different from 
the operating principles of present day MOS devices. 

We will now discuss these issues quantitatively by de- 
riving Eq.jSjfor gml I and evaluating its electrostatic and 
conformational components for specific illustrative exam- 
ples, with detailed numerical calculations for support. 

Transconductance equation. A source-drain bias Vd 
splits the contact chemical potentials /ii,2 by the applied 
voltage. The molecular I-V is obtained from Landauer 
theory [IS] , which can be recast in the less familiar form: 



/ = 



2e 



dE ( f{E) )■ 



f{E) = T{E)®Ft{E). (3) 

Thermal broadening appears jl4| as a convolution ® 
of the transmission T{E) at zero lead temperature with 
the thermal broadening function Ft{E) = —df{E)/dE, 
f{E) = l/[exp{E/3) + 1]. (...) gives a thermal average 
over various molecular configurations {x} with probabili- 
ties set by the conformational potential energies U{{x}): 



(TiE)) 



J2iT{E;Xi)wi 



w,=exp[-l3Uix,)], (4) 



Xi denoting the coordinate of the ith configuration. From 
Eq.|3 / — (2e//i)(/ii— /i2)(To), Tq being the awera(7e value 
oiT{E) over the energy range /ii < E < fi2. From Eq.Ql 
straightforward algebra gives 5m = ffm + 5™"^; where 
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(5a) 
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The transconductance has two additive contributions: an 
electrostatic term describing how the gate controls the 
channel charge, and a conformational term describing 
how the gate controls the transmission by deforming the 
molecule We now discuss these terms one by one. 

Electrostatic control. Electrostatically the gate modu- 
lates the self-consistent potential, effectively moving the 



molecular energy levels relative to the contacts. A com- 
plete treatment of this effect requires a self-consistent 
solution of the Schrodinger-Poisson equations. However, 
one can get a qualitative description by assuming that 
a gate bias rigidly shifts the energy levels and hence the 
transmission: T{E,Vg) « T{E — ^eVg). The parameter 
7 tells us the average molecular potential Vm in response 
to a change in gate bias: 7 = dVm/dVg. Using Ea.lKal 
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(6) 



7 ~ Cg/{Cg + Cs + Cd + Cq), where Cg,s,d represent 
the effective capacitances coupling the molecule to the 
gate, source and drain respectively, while Cq is the quan- 
tum capacitance proportional to the molecular density of 
states 17]. One could view Cg,s,d as a representation 
of Poisson's equation, and Cq as a linearized represen- 
tation of Schrodinger's equation. Under OFF conditions 
Cq is negligible and 7 can be made close to one in a well- 
designed FET by making iox ^ L, so that Cg ^ Cs.d- 
However, as we noted earlier, this becomes increasingly 
difhcult as we try to engineer ultrasmall devices that are 
only tens of atoms long [L ^ 1 — 5 nm) . The electrostatic 
parameter 7 is then severely reduced, motivating us to 
look for non-electrostatic mechanisms for control. Note 
that the electrostatic restriction on 7 is fundamental, and 
may not be easy to handle without using large dielectric 
constant insulators with a small effective thickness. 

Eq. El also elucidates the role of a sharp transition in 
To as a function of E in realizing good electrostatic con- 
trol. A metallic conductor with constant Tq{E) cannot 
be used to build an electrostatic switch. Semiconductors 
by contrast have a band-edge where To{E) drops to zero 
over a very small energy range AE, so that from Eq. |S1 



g'^^/I « ej/AE. 



(7) 



Since Tq{E) represents a convolution oiT{E) and Ft{E), 
the minimum value of AE equals that for Ft{E), which 
is ^ 1//3. This leads to the well-known upper limit of 
e/3 for gf,l/I, corresponding to 7 = 1 and AE = 1/(3. 
An important point to note is that increased tunneling 
in nanodevices can make T{E) non-zero below bandgap, 
with AE ^ 1//3. This condition is made worse by the 
use of metallic (rather than semiconducting) source and 
drain regions which is common for molecular devices P, 
0]. In view of recent progress in growing molecules on 
silicon surfaces , using doped semiconducting contacts 
seems a realizable and highly desirable ideal. 

Conformational control. The conformational mecha- 
nism (Eq. I5b|l operates by changing the relative ener- 
gies U for different molecular configurations {xi}. As ex- 
pected, we get zero transconductance if To is independent 
of the configuration {xi}, or its variation is uncorrelated 
with dU/dVg. But if Tq and dU/dVg are negatively cor- 
related such that configurations with larger transmission 



3 



To have their energies U reduced by the gate voltage Vg 
(and hence made more likely) then the current increases, 
with a positive 5™"^//. To estimate g^^^ /I, we write: 

U{e; Vd-, Vg) = Uo{e) - {iiVg/t,^) sine - {y.VjL) cose, 

. (8) 

where e is the angle between the molecular dipole fl and 
the source-drain field Vd/L. Aligning the dipole along 
the latter eliminates its torque [l^l, giving the gate field 
Vg /tox an obvious superiority. Using Eq. |H1 in Eq. I5bl 
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(9) 



The quantity within parantheses has a maximum value 
of one, indicating that the maximum conformational 
transconductance is given by /x/3/fox as stated earlier (see 
Eq. If /X is comparable to etox, then it is indeed pos- 
sible for this mechanism to provide respectable levels of 
control. For a 10 A oxide, n/etox 0.15 for an aro- 
matic molecule with one redox NO2 sidegroup per ben- 
zene ring. Although it is hard to squeeze in any more 
dipolar subgroups per ring, one could incorporate large 
effective dipoles external to the molecule, using a strong 
piezoelectric gate- molecular coupling, for instance In 
any case, a conformational transconductance that is, say, 
a tenth of e/3 may still be useful in view of the difficulties 
with the electrostatic mechanism discussed earlier. 

Significantly, the upper limit on the conformational 
transconductance (Eq.jSJ is unaffected by the 'softness' of 
the modes described by Uo{e) in Eq. |S1 For the bending 
and twisting configurations in Figsn^b,c), we write: 



Uoie) = u^''^" {e - eof /2, 

Uo{e) ^ [/**''^*(l-cos26l)/2. 



(10) 



Twisting a molecular bond is much easier than bending 
it (i/gt'^'id - 2 eV 0, f/o*^'"* ~ 0.03 eV |2l|); however, 
twisting also makes thermal averaging more significant, 
so that no advantage is gained as far as gm is concerned. 

The quantity inside parentheses in Eq. |51 can ideally 
have a maximum value of one if Tq is very sharply vary- 
ing. In practice it is smaller, depending on the variation 
of To with e. As an example, we model a self-assembled 
PDT monolayer (Fig.^) chemisorbed on a Au(lll) sub- 
strate, varying the tilt angle from default (^o ~ 20°) 
to upright position and increasing the sulphur-gold cou- 
pling exponentially. Although accidental symmetries at 
specific angles of the highly directional orbital wavefunc- 
tions quench the overlap , such orientational ef- 

fects wash out on averaging over various possible posi- 
tions of the sulphur and gold surface atoms in the upright 
configuration. The dominant angular dependence in To 
then arises from the overlap between sulphur and gold ra- 
dial wavefunctions: To(0) oc exp [—2ZL{\ — cos0)/naQ], 
{Z = 3: screened nuclear charge, n: principal quantum 



number of sulphur, oq = 0.529A). This gives a corre- 
lation in Eq. 1^1 of ^ 0.1 at room temperature (/Lt « 8 
Debye,K)/tox = IV /nm). This variation can be made 
to approach a delta-function by increasing T, giving a 
stronger correlation. Similarly for the rotational configu- 
ration in Fig. ^ we find Tq (0) cx cos'' e ji^ , leading to a 
configurational average ~ 0.3. The correlation depends 
on the variation To(0), but also on the weighting factors 
which in turn depend on the bond stiffness, temperature, 
gate field and dipole moment. 




FIG. 2: The I-V for PDT at a given angle saturates on cross- 
ing the highest occupied molecular orbital (HOMO) level, in- 
creasing with decreasing tilt angle. The electrostatic gate in- 
fluence (inset) causes a shift in the conductance peak towards 
lower voltages for p-type (HOMO) conduction. 

Numerical Results. We now present detailed calcula- 
tions of the I-V for the structures in Figs. ^QJ;. The 
transmission at a given angle is calculated using the 
nonequilibrium Green's function (NEGF) formalism |2^ 
with an extended Hiickel Hamiltonian H . We solve 

G{E) = (£;5-i/-C/(p)-I]i-S2)-\ 

/oo 
(/iGTiGt + /2Gr2Gt)rf£;/27r, 
-00 



T ^ trace (FiGFaGt) 



(11) 



self-consistently {S: overlap matrix, /i,2 = f{E — ^1^2))- 
An ideal Au(lll) surface geometry j24] is used to calcu- 
late the contact self-energies Ei^2 and level broadenings 
ri^2- The self-consistent potential U{p) is obtained by 
solving Laplace's equation for the device geometry with 
a 10 A oxide and a Hubbard-type electron-electron term. 

Fig.|21shows the zero-temperature I-V for a PDT relay 
for various tilt angles. The current tends to saturate on 
crossing the HOMO level, and increases with decreasing 
tilt (increasing Vg). In addition, there is a conventional 
electrostatic gate infiuence ji9j. For a p-type (HOMO) 
conduction, a negative gate bias raises the molecular lev- 
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els relative to the contacts, producing a lateral shift in the 
conductance (G = dl/dVd) towards lower \Vd\ (inset). 

Fig. |31 shows the three-terminal I-V of the Tour-Reed 
molecule |25| containing a nitroamine redox group at 
(a) K and (b) 300 K [2^. Ep is assumed to lie in 
the HOMO tail, giving an ohmic I-V at low V^. The 
impedance is compromised by the poor electrostatics due 
to the oxide thickness and MIGS, as discussed earlier. 
The low-temperature gm/I is impressive due to the low 
torsion constant and the large dipole, a mere 400 mV 
gate bias reducing the current two hundred times. At 
room temperature, however, the molecule samples a wide 
range of angles, reducing gm/I substantially. The gate- 
modulation continues to be observable, but the current is 
a lot harder to switch off. A lower torsion constant gives 

bpttpr cratp rnntrnl but inrrpasps 1:hprma.l pffprts as wpH. 




FIG. 3: FV of the Tour- Reed molecule )25[l for various gate 
voltages at (a) K and (b) 300 K, Ef lying in the tail of 
the HOMO level. At 300 K there is some gate modulation, 
although it is hard to completely switch the current off. 



The performance of a transistor also depends on its 
operating speed. While complicated isomeric cis-trans 
rotations tend to be slow, bond rotations are much faster 
(~ 10-100 GHz), and need to be damped out without 
consuming too much power. Image forces and Van der 
Waals interactions with electrodes jS^I as well as steric 
and hydrogen bonding interactions between molecules Q 
in a mixed monolayer could damp out such oscillations. 

We have shown that conformational transitions can aid 
electrostatic gate control significantly if we could engi- 
neer a large molecular dipole along a suitable direction. 
Ways around such design restrictions require going be- 
yond Eq. IHl employing other kinds of electromechanical 
gate-molecular coupling (piezoelectric for example), or 
using non-traditional, bistable potentials with large im- 
pulsive gate voltages lil . Eq. gives us a general way 
to quantitively compare the transconductance of these 
various field-effect mechanisms in molecular transistors. 
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